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Abstract 

We have developed a numerical code to study the evolution of distorted, 
rotating black holes. This code is used to evolve a new family of black hole 
initial data sets corresponding to distorted "Kerr" holes with a wide range 
of rotation parameters, and distorted Schwarzschild black holes with odd- 
parity radiation. Rotating black holes with rotation parameters as high as 
a/m = 0.87 are evolved and analyzed in this paper. The evolutions are 
generally carried out to about t = lOOM, where M is the ADM mass. We 
have extracted both the even- and odd-parity gravitational waveforms, and 
find the quasinormal modes of the holes to be excited in all cases. We also 
track the apparent horizons of the black holes, and find them to be a useful 
tool for interpreting the numerical results. We are able to compute the masses 
of the black holes from the measurements of their apparent horizons, as well 
as the total energy radiated and find their sum to be in excellent agreement 
with the ADM mass. 

PACS numbers: 04.30.-hx, 95.30.Sf, 04.25.Dm 
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I. INTRODUCTION 



In this paper we present the first results for the evolution of initial data sets correspond- 
ing to axisymmetric, rotating vacuum black holes. We developed a numerical scheme and 
code based on an extension of earlier work at NCSA [0-^ on distorted, but non-rotating, 
axisymmetric black holes. In a companion paper, referred to henceforth as Paper I 0, we 
describe the numerical methods, gauge conditions and code tests developed for a code to 
evolve rotating black hole spacetimes. With this new code we have been able to evolve and 
study a new family of distorted rotating black hole data sets. The construction of these data 
sets was outlined in Paper I, and will be discussed in another paper in this series which 
we refer to as Paper III. In this paper we focus on analyzing the physics of the evolution for a 
number of different rotating black hole data sets, with both even- and odd-parity radiation. 
We also consider the evolution of non-rotating black holes that have been distorted by the 
presence of odd-parity gravitational waves. 

These rotating black holes can be highly distorted, allowing one to study their nonlinear 
dynamics. These studies will be useful not only in understanding rotating black hole space- 
times, but also in studying the late stages of black hole collisions with angular momentum, 
just after the holes have coalesced. At that time they will have formed a single, highly 
distorted rotating hole similar to the configurations studied here. 

In analyzing these rotating black hole evolutions, we have developed a series of tools that 
allow us to study gravitational waves, apparent horizons, and other quantities of interest. 
Using these tools we have extracted the waveforms for both the even- and odd-parity radia- 
tion emitted by a distorted rotating black hole, and we find that the quasinormal modes of 
the hole are excited. The extracted waveforms are also used to compute the energy carried 
away from the black hole via gravitational waves. We locate and study the geometry of the 
apparent horizons in these spacetimes, and find that their geometric structure can be used 
as a key tool in analyzing the evolution. We show how to extract information about the 
oscillation frequency, the mass, and the angular momentum of the black hole. 

The organization of this paper is as follows: In section we define the variables used 
in our simulations and describe the variety of initial data sets we evolve and analyze in this 
paper. In section |T| we present an analysis of "near zone" results from evolving various 
initial data sets, including an analysis of the physics that can be extracted from apparent 
horizons. Section |^ contains a discussion of our extraction methods to compute the even 
and odd-parity radiation in the spacetime, while section ^ contains a discussion of the results 
of evolving the data sets. Finally, in section [V^ we conclude and discuss future directions 
this research will take. 



II. THE CODE AND INITIAL DATA 

The code we developed to evolve rotating black hole data sets has been discussed in 
detail in Ref. We will not discuss any details of the numerical methods used to evolve 
the spacetimes here, but for context we define our notation for the primary variables used 
in our code. 
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A. Definition of Variables 



We build on earlier work of Ref. [|T],^ on non-rotating black holes in defining the variables 
used in our code. Additional metric and extrinsic curvature variables must be introduced 
to allow for the odd-parity modes present in this new system. Previously only even-parity 
modes were present. We define the variables used in our evolutions as follows: 



and 



IrjTi IrjO Irjcp 

lij = lr,e lee le^ 



I A C Esm^e 
C B F sin e 
V^sin^^ Fsin^ Dsin^^ 



(la) 



Ha He He sin^ 9 
He Hb Hp sin 9 
\ He 3111^9 Hp sin 9 HDsin^9 



(lb) 



In these expressions rj is a. logarithmic radial coordinate, and {9,(f)) are the usual angular coor- 
dinates. The relation between t] and the standard radial coordinates used for Schwarzschild 
and Kerr black holes is discussed in Paper I. As in Ref. [Q, the conformal factor is de- 
termined on the initial slice and held fixed in time afterwards. The introduction of into 
the extrinsic curvature variables simplifies the evolution equations somewhat. The various 
factors of sin^ are included in the definitions to explicitly account for their behavior near 
the axis of symmetry and the equator, as discussed in Paper I. 



B. Initial Data 

We constructed a new class of initial data sets corresponding to distorted rotating black 
holes and odd-parity distorted non-rotating black holes. The construction of these data 
sets is described in Paper I and in more detail in Paper III. Here we briefly describe the 
different classes of data sets and describe our parameter choices for evaluation and analysis. 
There are three basic families of black hole data sets that are evolved and analyzed in this 
paper. The first series of data sets, labeled in Table | as runs kl and k2, correspond to 
Kerr black holes with rotation parameters a/m = .48 and .68 respectively. As described in 
Paper I, these data sets have been transformed from Boyer-Lindquist coordinates into our 
logarithmic coordinates and then evolved. 

The second family of initial data sets that we consider are related to the rotating black 
hole initial data of Bowen and York 0|, but we have added a "Brill wave" to distort them 
with a construction similar to that of Ref. 0. As described in Paper I, the initial three- 
metric 'Jab is chosen to be 

where the "Brill wave" function q is given by 

q = sin" 9 qc, (3a) 

qG = Qo{e-'^ + e-'-), (3b) 

s± = {r] ±r]of /a"^. (3c) 
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The extrinsic curvature is computed following Bowen and York, with a free parameter J that 
specifies the angular momentum of the hole. For more details of how the constraints are 
solved, please refer to Paper I. Here we simply point out that the Bowen and York solution 
to the momentum constraint is: 

He = 3^V, (4a) 
Hf = 0. (4b) 

We note that for J = 0, the data sets reduce to the "Brill wave plus black hole" distorted 
black hole data sets described in Ref . [|1| , and for the Brill wave parameter g = the Bowen 
and York solutions result. 

The initial data sets described in Table |1| whose name begins with an r {rO through r5) 
describe a sequence of distorted rotating holes of increasing rotation. The run rO has no 
angular momentum, while run r5 is a pure Bowen and York black hole without a Brill wave 



added, but with a significant amount of rotation. As we will see in section |III A 3| below, it 
is rotating so fast {a/m = 0.87) that its horizon cannot be completely embedded in three 
dimensional Euclidean space, yet we are able to evolve it accurately. 

Finally we have considered a third family of distorted black hole data sets that correspond 
to odd-parity radiation superimposed on non-rotating black holes. In this case the three- 
metric Jab is described by the same Brill wave parameters in Eq.(^ above, but the extrinsic 
curvature is taken to be 



He = qc [{n + 1) - (2 + n') sin^ 6) sin"'-^ 6 (5a) 
He = -^^ d^QG cos e sin"'-i 6. (5b) 



The parameter n' is used to describe an "odd-parity" distortion. It must be odd, and have 
a value of at least 3. This family of data sets does not describe rotating black holes, as 
discussed in Paper I and in more detail in Paper III. The two data sets in Table || whose 
name begins with o (ol and o2) represent a family of odd-parity distortions of non-rotating 
holes. As we will see in section |III A 3 below, the angular index n' is related to the dominant 
i-mode in the distortion of these holes. Both an n' = 3 and an n' = 5 data set are given 
because the former emitted virtually all of its radiation in the £ = 3 mode and we were also 
interested in examining significant i = 5 radiation from an odd-parity distortion. 

We summarize the dimensionless angular momentum parameter {a/m) of these space- 
times in Table ^ where three different measurements are tabulated. The meanings of the 
columns are as follows: {a/m)min = J/^adm^ would be the rotation parameter of the 
system if none of the radiation escapes to infinity. {a/m)max = J/^ahi where Mah is the 
mass of the apparent horizon defined by Eq.([10|) below, and would be the rotation parameter 
of the final system if the area of the black hole did not increase at all during the evolution. 
These two measurements are made on the initial data sets without evolving them. Finally 
{a/m)cr is the value of {a/m) obtained by measuring the shape of the apparent horizon 
surface during the evolution as discussed below and in Ref. 0. 
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III. NEAR ZONE PHYSICS 



In this section, we discuss details of the "near zone" evolution of several black holes, 
focusing on the apparent horizon. We present an apparent horizon finder, and discuss the 
dynamics of the apparent horizon as a tool to study the physics of the near zone. Although 
extremely important to the numerical evolution, the behavior of the metric functions them- 
selves generally does not illuminate the physics of the system, and is not discussed here. We 
refer the interested reader to Paper I for discussion of metric evolution. 



A. Horizons 



In this section we examine the properties of the black hole apparent horizons in the 
spacetimes evolved here. The apparent horizon is defined as the outermost trapped surface. 
Because a knowledge of the future is required to identify the event horizon, the concept of 
an apparent horizon is useful in numerical relativity. The location of the apparent horizon 
depends only on information known within a given time-slice, is always inside the event 
horizon [ffO] and it coincides with the event horizon in a stationary spacetime. As shown 



in Ref. ||Tl|, where studies of perturbed event and apparent horizons were made, apparent 
horizons can closely approximate the event horizons even in dynamic spacetimes, and share 
many dynamical features with them. In this paper we focus on the apparent horizon, but 
will treat the event horizon of highly distorted black holes in a future work. 

Mathematically the apparent horizon condition is expressed by requiring the expansion 



of all outgoing null normals to the surface to vanish. This condition can be written as [|12 



Das" + KabS^s' -trK = (6) 

where s° is the outward pointing normal vector to the surface in the three dimensional 
subspace. If such a surface exists, it is said to be a marginally trapped surface. In some 
cases multiple surfaces satisfying this condition can be found, and the outermost one is then 
defined to be the apparent horizon. The apparent horizon finder searches for a solution to 
this equation by making an initial guess, and then iteratively solving the above nonlinear 
equation via the Newton- Raphson method until we obtain the desired accuracy 

The initial guess that we use is an t] =constant surface where the value of 1] used is the 
outermost value which satisfies the condition 

on the equator [9 = 7r/2). This may be interpreted as the condition that the light rays 
confined to the equatorial plane have null-expansion. The radius thus determined provides 
us with a good estimate of the position of the apparent horizon. This technique proved 
valuable because it made it easier to follow the horizon when it jumped outward from the 
throat at about t ~ lOM. This particular complication is a result of using an antisymmetric 
lapse that vanishes on the throat of the black hole, preventing evolution there (except 
through the shift terms). The initial apparent horizon is usually located on the throat (the 
throat is guaranteed to be a minimal surface, but not necessarily an apparent horizon). With 
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a lapse that vanishes on the throat, the horizon often remains there until some evolution has 
occurred, at which point it may jump out discontinuously. This technique does not require 
the horizon position from previous time steps to be recorded, allowing us to find the horizon 
as infrequently as we wish, or to find it again if it becomes lost. This trick has been quite 
successful and uses up little CPU time. Once we have located the horizon with this method, 
it can be analyzed to provide physical information as we describe in the following sections. 



1. Masses 

In addition to the ADM mass, which was calculated on the initial slice, we are interested 
in obtaining a dynamic mass during the evolution. The measure we will use for this is 
an "areal mass" based on the irreducible mass (as given by the apparent horizon) and the 
angular momentum of the black hole. 

The horizon finder returns the location, rjiO) of the apparent horizon on a given slice. 
In order to study the geometric properties of the horizon, from which important physical 
information can be obtained, it is useful to determine the two dimensional sub-metric induced 
on the surface by defining a new coordinate ^ by the equation i] = t]{6) ^. We now have 

df = (^Ad^^ +^B + A^^ (^^(^)) j d9^ + Dd(f + 2Fd<j)de^ . (8) 

The sub- metric that we are interested in may be found by letting ^ be equal to 1. The 
surface area of the horizon is, therefore, 



A = 2tx siYiO^^dd 







\ 



For a Schwarzschild black hole the mass is given by Mj^ below, and for a Kerr black hole 



the mass is determined by both the surface area and angular momentum of the hole [jlS 
according to the relationships 



and 



Although our spacetimes are not stationary and hence the apparent horizon does not coincide 
exactly with the event horizon, these are quite useful quantities to take for the dynamical 
mass. The above formulations describe the minimum mass of the final Kerr black hole that 



is possible at the end of the evolution However, as shown in Refs. |]Tl] , p!5| the event and 
apparent horizons are often quite close, even in dynamic spacetimes, and the mass described 
in Eqs. ( pUD should often be a good approximation to the event horizon based mass. The 
masses computed in this way will be used extensively in analyzing the evolutions presented 
in section below. 
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2. Oscillations 



Once we have the 2-metric on the horizon's surface we can measure the horizon geometry. 
An example of such a measurement is Cr, the ratio of polar circumference (Cp) to the 
equatorial circumference [Cg) of the apparent horizon. The horizon of a stationary black 
hole has a characteristic shape that depends upon its rotation parameter a/m [|T6|. For a 
stationary or a dynamic black hole, this shape can be used to obtain the rotation parameter, 
as discussed in Ref. |0. There it was shown that the horizon oscillates about this equilibrium 
stationary shape with the quasinormal mode frequencies of the black hole. We can plot the 
shape parameter Cr as a function of time to measure these horizon dynamics. Here we 
extend this technique and apply it to new spacetimes not considered in Ref. 0. 

To obtain the polar and equatorial circumferences we first transform the two-metric so 
that it becomes diagonal. This is accomplished by introducing a new coordinate, x- The 
metric on the two surface defined by the apparent horizon now becomes: 

F 

dx = d(f)- . d9 (11a) 
U sm (J 

F' = (11c) 
D' = m*D. (lid) 



The polar circumference and the equatorial circumference may now be defined as: 

Cp = jde^, (12a) 

and 

Ce = <fd(j)y/D'. (12b) 



It was shown in previous work that the i = 2 and i = 4 quasinormal mode (QNM) 
frequencies can be seen in the function Cr, which can also be analyzed to obtain the rotational 
speed [a/m) of the spacetime. This is accomplished by first fitting Cr to the fundamental 
£ = 2 QNM, the first overtone i = 2 QNM, and a constant offset from unity. We assume 
that the oscillations are dominated by the i = 2 quasinormal ringing radiation going down 
the horizon. As discussed in Ref. the ringing radiation is generated at the "peak of the 
potential barrier" located outside the horizon. Part of this signal propagates away from the 
hole and part goes down the hole, causing the horizon to oscillate as the signal crosses it. 
The offset is used in the fit because a rotating black hole will not be spherical (a sphere 
would have Cr = 1) but should be oblate (in equilibrium), producing an "offset" in Cr- The 
fitted value of the constant offset can then be used to obtain the value of a/m for the surface, 
since for a Kerr black hole the oblateness is a unique function of the rotation parameter. 
(We note that the following function is an approximation generally accurate to within 2.5% 
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for a/m as a function of Cr- a/m = y 1 — (—1.55 + 2.55 C^) .) In Table fT| we show the 
value of the rotation parameter a/m for the family of data sets labeled r0-r5 obtained with 
this method. 

Finally, as the horizon oscillations should be caused by radiation from all black hole 
modes excited during the evolution, not just the C. = 2 modes, one can analyze the function 
Cr for other modes as well. For example, the fit to the shape parameter in terms of the 
1 = 1 QNM expansion functions can be subtracted from the actual function C^, leaving a 
residual of higher modes. In Ref. we showed how this remaining signal was dominated 
by £ = 4 QNM's for some spacetimes. 

For the odd-parity distorted holes (oi and o2^ not shown), the oscillations in Cr are 
second order in the metric perturbation. As a result the spectrum of frequencies in the 
surface vibrations that result from the £ = 2 and £ = 3 QNM's will be much larger. Therefore 
the oscillations are more complicated than in the previous cases. Further details of the 
surface oscillations of these black holes will be given in section [lllB . 

In Fig. ^ we plot C^ as a function of time for simulation rJ^. This is a distorted Bowen 
and York hole rotating quite rapidly. The solid line shows the value of Cr extracted from 
the horizon. The dashed line shows the fit to the two lowest 1 = 1 QNM's, and the straight 
horizontal dashed line shows the constant offset that came from the fit. Note that the fit 
matches the oscillation quite well except at very late times, when we expect the code to 
be less accurate, particularly near the horizon where metric functions tend to develop steep 
gradients . This analysis shows that the black hole is oscillating about a Kerr black hole 
with a rotation parameter of a/m = 0.70, consistent with estimates made on the initial slice. 

In Fig. ^ shows a similar plot for the function Cr and its fit for run r5. Run r5 differs 
from runs rO-rJ^ in that r5 is a pure Bowen and York black hole. It is thus much more 
noticeably oblate, with Cr ~ 0.81, as opposed to Cr ~ 0.89 for the run labeled r^. In this 
case the analysis shows the rotation parameter to be a/m = 0.87. 

Table |T| gives the results of extracting a/m by this method in column 3 for runs r0-r5. 
This table shows that in all cases the value extracted from Cr during the run is very nearly 
the lower limit, J/M\j;)j^j. This means that in the runs under consideration in this paper 
nearly all the radiation is going into the hole, because the hole achieves nearly the maximum 
possible mass (or minimum possible rotation parameter). 



3. Emheddings 

A valuable feature of the horizon geometry is its three-space embedding. By embedding 
the horizon we can obtain a visual impression of its overall distortion in a coordinate inde- 
pendent way. The Kerr geometry is a sphere in Boyer-Lindquist coordinates, but it is not so 
when embedded in three-space. Instead it is an oblate spheroid that becomes increasingly 
flattened along its axis of rotation as its spin increases. If its rotation parameter a/m exceeds 
-\/3/2 it is impossible to completely embed the figure in fiat space, as shown by Smarr |]16 



The procedure for finding the embedding starts by defining a new coordinate z which is 
part of a flat 3-metric and is identified with the 2-metric of our surface: 

dz^ + dp^ + p^dcj)'^ = B'de^ + O'dcj)"^. (13a) 

Solving this expression for z gives 
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d9\ B 



{de 



'D' 



(13b) 



and 



P 



D. 



(13c) 



This equation is then integrated numerically to obtain the embedding functions z{9) and 
p{d). In Fig. ^ we show the embedding of the evolving horizon for run r^. At time t = 5 the 
horizon is still frozen and has its initial prolate shape with Cr = 1.228. At time t = lOM 
the horizon has become noticeably oblate, overshooting its equilibrium value, and by the 
late time of t = 50M it has settled down to very near its equilibrium Kerr value. The dark 
solid line shows the shape of the Kerr hole with the rotation parameter determined from 
the procedure above. (The area has been normalized to take account of the difficulty of 
resolving the peak in A that can cause the area to grow spuriously, as discussed below and 
in Ref. §.) 

For a Kerr black hole, when the rotation parameter a/m becomes large enough and the 
integrand of Eq. (|13b| ) becomes imaginary, the embedding ceases to exist. However, this 
generally happens only over a small portion of the horizon surface, near the symmetry axis. 
This embedding failure will begin to occur when a/m = V3/2 ^ 0.87M. The angle 9 
(measured from the z-axis) at which the embedding should cease to exist can be found. The 
equation for this angle is 

2n 5 + 36 4 2a;V3 

cos & = — jtt:; 77^, (14a) 

3(6-1) 3(l-6)'/'a;i/3 3(1-6)'/=^ ^ ^ 

where 

X = 17 + 106 + 3y^37im26T462), (14b) 

and 



{a/mf. (14c) 



Note that for the region 6 ^ Att the embedding never disappears. 

This effect can actually be used to measure the rotation parameter for rapidly rotating 
holes. For a dynamic rotating black hole, as shown above the horizon geometry will oscillate 
about the stationary Kerr shape. For an extremely fast rotator or highly distorted black hole 
the oscillations may distort the horizon so much that the embedding fails at times during 
the evolution (or in extreme cases the embedding may never exist). In Fig. ^ we plot the 
value of 6 for which the embedding fails as a function of time for the case r5. This is a pure 
Bowen and York hole with J = 15, corresponding to a/m = 0.87. 

The run labeled r5 is the fastest rotating black hole we have simulated to date, and it 
pushes the limits of our code, but it nevertheless produces a good result, in agreement with 
theory. The straight horizontal line represents the embedding limit of the analytic Kerr 
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solution with the given value of a/m = 0.87, which was the rotation parameter that was 
extracted from the oscillation about Cr as discussed above. Note that this black hole horizon 
is never globally embedded in three dimensional Euclidean space, although the horizon 
always exists. Note also that by measuring the horizon geometry, we can determine the 
rotation parameter at least two different ways: by measuring Cr as in Fig. |1| (see discussion 
above), and by measuring the angle at which the embedding ceases to exist (although this 
only works for holes which have a/m > ^/3/2). Furthermore, these measurements agree 
with estimates of the rotation parameter determined from the initial data alone. 



B. Horizon History Diagrams 

As we have shown in Ref. it is useful to construct a "horizon history" embedding 
diagram, that shows the evolution of the surface in time. In order to bring out details of 
the local curvature of this surface, we also compute the Gaussian curvature k of the horizon 
surface and map it to a colormap or greyscale. For details of this construction we refer the 
reader to Ref [Q, but for reference we present the formula for Kerr here: 

8 (-1 + 56-3 cos(20) +36 cos(20)) 
^ ~ (1 + 6) m2 (3 + 6 + cos(2 6) - b cos(2 e)f ' ^ ^ 

where b = y'l — /w? as before. The mapping of the Gaussian curvature onto the surface 
helps to bring out small deviations in the local curvature of the surface that would not be 
apparent in the embedding diagram itself. 

In Fig. 1^ we show a horizon history embedding diagram for the case oi, which is an 
odd-parity distorted non-rotating black hole. As we discuss in detail in section 0, this data 
set has predominantly £ = 3 radiation (99% of the total energy radiated is carried in this 
mode). As in the cases shown in Ref. each £-mode pattern is qualitatively different from 
the ones studied previously. In Fig. ^ we show a similar diagram for the run labeled o2, 
which has a significant £ = 5 component. Again, the pattern is qualitatively different from 
the £ = 3 (or £ = 2 or £ = 4) pattern, and has not been seen before. 

These distinctive patterns can be understood by analyzing the expression for the gaussian 
curvature. The metric variable F enters the Gaussian curvature as a second order term (to 
first order it vanishes identically). This function carries the odd-parity radiation in our 
gauge, as discussed below. Therefore if F oscillates predominantly at the £ = 3 normal 
mode frequency denoted by uj^=^, then the period of oscillation one will see in the Gaussian 
curvature plot will not be a;£=3 due to this nonlinearity. Instead one will see 2uj^=^, so the 
pattern cycle will repeat twice as quickly as it would have if F had entered the Gaussian 
curvature to linear order. This is what one sees in Fig. ^. The pattern repeats every 
if: 5.2AM instead of every t ~ 10.48M, which is the period of the £ = 3 fundamental mode. 

In Fig. P the i = 5 modes are present to a much stronger degree, so we see a mixture of 
the i = 3, and £ = 5 modes. This creates a diamond pattern reminiscent of the diamond 
pattern seen on the surfaces of even-parity black hole horizons generated with an n = 4 
perturbation, as described in Ref In general, however, this pattern is more complicated 
than the n = 4 pattern, as it should include four frequences: 2uJi=3, u!£=3 ± co'£=5, and 2u!i=5. 
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IV. FAR ZONE PHYSICS 



A. Wave Mode Extraction 

In previous sections we have analyzed properties of the near zone features of these new 
black hole data sets, focusing on geometric measures of the horizon. These effects provide 
important probes of the physics of the near zone, and can also be used to determine important 
properties of evolving black holes that may not be known, but for the most part these effects 
are not measurable by gravitational wave experiments such as LIGO. In this section we turn 
to analysis of the physics of the far zone, where important features such as gravitational 
waves can actually be observed. 

One of the principal features of a dynamical black hole spacetime which has been possible 
to study analytically is linearized gravitational waves. Therefore, it is useful for us to measure 
them in highly distorted black hole spacetimes to determine information that linear theory 
cannot predict, such as the final mass of a black hole in a perturbed spacetime, or the 
waveform emitted during the nonlinear generation of gravitational waves, and to provide 
several useful checks of our code. Although the processes that generate these waves may be 
highly nonlinear, far from the hole the waves may be treated as linear perturbations on a 
fixed background. For all the simulations presented in this paper we extracted the radiation 
at a distance r = 15M. 

The radiation energy in a black hole spacetime is described by certain gauge-independent 
variables constructed from the metric and its derivatives. The technique was originally 
developed by Abrahams and Evans and applied to black hole spacetimes in Ref. 0. 



There are two classes of gauge-invariant radiation quantities representing the two degrees of 
freedom of the field, even-parity and odd-parity. Both are present in our spacetimes and we 
describe them below. 



1. Even-parity 

We used the same wave extraction routines for the even-parity wave forms as used in 
Ref 0]. The routines were developed to extract waveforms from perturbed Schwarzschild 
black hole spacetimes that oscillate about a spherical background. Because the distorted 
Kerr metric does not settle down to something that is spherical, but rather it settles down 
to something oblate, we should expect to see this effect in the extraction process, resulting 
in an offset in the i = 2 and £ = 4 waveforms. Thus, the wave does not oscillate about 
zero, but instead will be offset from zero depending on both the rotation parameter and the 
radius at which the wave is extracted. 

However, the level of this offset is small and its variation rather too sensitive so that we 
have not been able to use it to measure the rotation parameter reliably. As shown below, 
we do see this effect in waveforms extracted from distorted rotating black holes, and the 
offset value is related to the rotation parameter as expected. The offset we obtain from the 
waveforms is of the right general magnitude and sign for a given rotation parameter. In 
principal one can account for the fact that the system is a perturbed Kerr black hole, but 
to date we have not carried out this analysis. However, as we show below, the extraction 
assuming a spherical Schwarzschild background can be quite useful without modification. 
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The even-parity wave extraction is given here. Note that this extraction formula as- 
sumes a Schwarzschild background. For each £-mode we can extract independent radiation 
waveforms. The gauge invariant, even-parity wave function is given by 



\ 



2{£ -!)(£ + 2) {4:rS^k2 + i{i + l)rh] 



and 



+ 11 



A = 1) -2 + 



A 



6M 



(16a) 



(16b) 



where the Moncrief functions are given by 

dG 



ki = K + Sr 
ko = — 



dr ' 
1 d 



2S 2S^2 dr 



rS-^K 



with 



5 = 1- 



2M 



(16c) 
(16d) 



(16e) 
(16f) 



and finally the Regge- Wheeler |2^ perturbation functions are defined in terms of the 3-metric 
via 



H2 



5 / m^AYtosinede, 



A^ Jo 
1 



2 Jo 



sinOdd, 



K 



i?2 70 

£{£+!] 



+ ii + 2) (i-l) 

+ 4 (B + D) sin e Ymde, 

K"^ Jo 



1 

2 Jo 



sin^ci^. 



(16g) 
(16h) 
(16i) 

(16j) 
(16k) 



With the normalization of ip given above one can show that the total radiated energy in 
each £-mode is given by 

E = ^jdt id.ijf . (17) 
Complete details of this extraction procedure are provided in Ref. 0|. 
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2. Odd-parity 



In terms of the Regge and Wheeler formahsm PO] there are various ways to construct the 



odd-parity gauge-invariant variable that measures the radiation in the system ||19[. These 
measures are not linearly independent, and each is constructed from two of the variables 
E, F, and (3'^, which correspond to Regge and Wheeler's variables hi, h2, and ho after 
appropriate angular integrals have been performed. The gauge-dependent Regge- Wheeler 
variables can be extracted from the metric as 



ho 
hi 



(i+lV r 



£0} 



[18a) 
;i8b) 



and 



(£ + 2) ! 
2^ (£-2) 



- j dn ^^F {de - cot 9) deYa 



;i8c) 



From these perturbation functions one can construct the following gauge-invariant quanti- 



ties, two of which are given in Ref |21 



i^02 = ho + -dth2, 
ijjoi = r'^dr ( ^ I - 



dthi, 



(19a) 
(19b) 



and 



= M + ir=9, ^1 . 



(19c) 



In our code we have used our gauge freedom to eliminate the metric function E (hi) 
and the method we use to extract the radiation is effectively Eq.( |19c|) . It is also possible 
to measure the odd-parity radiation through Eq.(r^) as we show below. Note that for the 



odd-parity case, there is no "non-spherical odd-parity" part in a Kerr black hole (in iIji2) so 
there is no offset expected nor observed for these waveforms. 

The particular energy integrals for the odd parity modes are given here. Note that these 
expressions are specific to our gauge and a flat space background metric. We first normalize 
our gauge invariant odd-parity wave function, based on Eq. (|19c|) as follows: 

^°^^3 = £ de ^VIOStt cos 6 sin^ OdnF (20a) 

and 
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r rf^— VllSSvr (5 cos ^ + 3 cos (3^) ) sin^ Odr^F. (20b) 
JO 16 



For general ^ the energy expression is this 

^odd = r 27r. 2^^f^de d^F sin 9 {de - cot 9) deYto- (21) 
JO \ (i; + 2j! 

In principle one can extract arbitrarily high £-modes. In practice we have only examined 
the £ = 3 and £ = 5 modes to date. With these normalizations, the energy radiated by each 
£-mode is given by 

^=3^ (22) 



as in Ref. 21 . 



B. Measurement of Frequencies 

To analyze quantitatively the frequency spectrum of our radiation we used the following 
Fourier technique. Because the quasinormal modes of black holes are all damped the Fourier 
spectrum of their frequencies will be spread out. To counter this effect, before analyzing a 
waveform we premultiply it by a factor exp(At), where A is chosen to approximately cancel 



the effects of the damping (see Ref p2| for a discussion of such techniques). As it happens, 
the damping time for the fundamental mode is given approximately by A ~ .09/M, relatively 
independent of the £-value, so this technique is helpful for all modes. 

We next decompose the frequency spectrum using the procedure described below. The 
technique used here to analyze the waveforms is preferable to an FFT, since an FFT may 
provide poor frequency resolution for the wavelengths under study unless the sampling 
interval (AT) is quite long. For example, to distinguish between the peaks in ^ = 2 and 
£ = 3 QNM's would require a "sampling period" greater than AT 60M of oscillations. 
The Fourier technique we present below, combined with the premultiplication technique, 
can distinguish between them with only about AT ^ 17M of oscillations. We approximate 
the Fourier transform function and its inverse by 

fiuj) = -= / f{t)e~^-'dt ^ f^^^i^)^t (23a) 

V ZTT J-oo j 

fit) = -= / f{u)e^-'du ^ /^^(^ - i^At, (23b) 



where we have defined 



and 



6{t-U) = ^e~^'-'^^"/'^\ (23c) 



u) = e-*-*>-'-V4_ (23d) 
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In the above equations and f(t) is the Regge- Wheeler or Zerilh function, {fi,ti) are the set 
of points produced by the numerical code (we should find that /, ~ f(ti)), and t is the time. 
The value of a was set to At. 

Note that the presence of the first overtone of a quasinormal mode frequency cannot be 
as easily detected with this technique (since the damping is so much greater). It manifests as 
a small shift in the position of the peak of the fundamental mode. The position of the peak 
will also be affected by the artificial growth in the mass that occurs at late times, which is 
manifested in the waveform by lengthening of the wavelength of the radiation as discussed 
in Ref. 0. This artificial lengthening will effectively add a mixture of other frequencies to 
the spectrum, further broadening the peak. 



V. RESULTS AND DISCUSSION 

In this section we discuss results from evolutions of the initial data sets listed in Ta- 
bles |I| and ^ combining many of the analysis techniques discussed in the previous sections. 
Together these tools provide a thorough and remarkably consistent physical picture of the 
evolution of these black hole spacetimes. 



A. Odd-parity Distorted Schwarzschild Black Holes 

In this section we give results for evolution of the new class of data sets we constructed 
and discussed briefly in Ref. 0. These evolutions correspond to runs labeled ol and o2. 
These data sets do not possess angular momentum, but they do possess odd-parity radiation, 
in contrast to the "Brill wave plus black hole spacetimes" discussed previously in Ref. [|I|. 
Strictly speaking, they are odd-parity distorted Schwarzschild black holes, and not rotating 
holes. Nevertheless, the odd-parity distortions give rise to "rot at ion- like" features. We note 
that these black hole data sets also contain even-parity radiation, although of a much lower 
amplitude. 

The run labeled ol represents a spacetime with virtually only an £ = 3 distortion, as 
almost all the energy is radiated in that mode. Fig. ^ shows the i = 3 waveform extracted 
by the gauge-invariant waveform extraction method described in section [IV A2| above. As 



one can see, it is almost impossible to distinguish the fit from the data. As we show in 
Table III, 99.96% of the total energy radiated is carried by the £ = 3 mode. In Fig. | we 
show the i = 5 waveform extraction for the same run. It is interesting to note that even 
though it does not contribute significantly to the energy, and the signal amplitude is nearly 
four orders of magnitude smaller than the i = 3 signal, the i = 5 signal is still easily fit to 
the proper mode. 

For run ol we extracted the £ = 3 radiation mode from both (3^ and F. Both extractions 
are plotted in Fig. The solid line illustrates the extraction from the shift, the dotted 
line traces the value extracted from F. The waveform extracted from f3'^ was normalized so 
that it had the same lower bound as the waveform extracted from F. This shows that it 
is possible to perform the extraction from the shift, although we have only done so for this 
test. 
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Next we consider the run labeled o2, which has a larger £ = 5 signal. The i = 5 



extracted mode, shown in Fig. |10|, matches the fit to the quasinormal modes quite well, 
as in the previous case. (The i = 3 waveform, not shown here, matches its quasinormal 
mode fit with the same level of accuracy as in the previous case ol.) After measuring the 
i = 5 frequency from our extracted waveform via the Fourier transform technique discussed 
in section |IVB| above, we discovered that its quasinormal modes were not tabulated in 
the literature. We computed the real part of the frequency to be 1.00 ± .01. Subsequent 
calculations by Edward Leaver, based on black hole perturbation theory yield the result 
1.012 which agrees with our result to within 1%. In Fig. |TT| we compare the analytic i = 5 
frequency with the Fourier transform (as implemented by the technique described above) of 
our data for the run labeled o2. A dotted line is placed on the graph to show where the peak 
should be for a pure i = 5 wave. The secondary peaks do not represent real frequencies, 
rather they are an artifact of the extraction process. 

Finally, we turn to a calculation of the energies and masses in the system. The mass of 



the apparent horizon for the run labeled ol is plotted in Fig. 12. The solid line shows the 



instantaneous mass of the horizon, defined by Eq.(|TO|), and the long dashed line represents 
the total energy radiated away from the black hole during the evolution. All energies in 
this plot are normalized to the ADM mass, and therefore it can readily be seen that the 
sum of the final apparent horizon mass (as measured at t = 25M) and the energy emitted 
through radiation (as calculated by integrating the even- and odd-parity radiation functions 
for £ = 2,3,4 and 5 for the entire time of the run) add up to the ADM mass. The very 
slight dip in the mass at around 30M becomes less noticeable with increasing resolution. 

After about t = 35M the peak developing in the radial metric function is not adequately 
resolved. As the horizon is located near this peak, as discussed in Paper I, its area is 
not accurately computed after this time. However, the simulation can be continued until 
about t = lOOM and accurate waveforms can be extracted throughout. Furthermore, the 
calculations of the shape of the horizon continue to be accurate throughout the evolution, 
even though the area becomes inaccurate. 



B. Kerr Black Holes 

Although a Kerr black hole is stationary and not dynamical it does provide a useful test 
case for evolution. Furthermore, as it is the "equilibrium" black hole about which all other 
black hole spacetimes considered in the paper oscillate, it is useful to study its properties. 

Because there is no radiation in the system we do not expect the horizon to oscillate. 
In Paper I we used Kerr as a test case to show that our code was able to evolve the Kerr 
spacetime accurately by computing its angular momentum during the evolution. Here we 



apply two of the near zone measurements developed in section |T| for these spacetimes to 
show the numerical properties of these evolved black holes. 



In Fig. we show the mass Mah of an a/m = .48 black hole (labeled kl) as a function 
of time for three different resolutions. The horizon mass should be strictly conserved and 
equal to the ADM mass for a stationary black hole. The high resolution (300 x 30) mass is 
plotted as a solid line, the medium resolution (150 x 24) mass is plotted as a dotted line, and 
the low resolution (75 x 12) mass is plotted as a dashed line. The mass is nearly constant 
until the large peaks discussed in Paper I develop and cannot be adequately resolved. By 
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t = 40M, the error in the mass is still less than .1% at all but the lowest resolution and by 
t = 60M the error is still only .5%. Although the apparent horizon mass does not exactly 
equal the ADM mass, it is within the error expected given the level of angular resolution of 
the horizon. The jaggedness of the lower resolution lines is a numerical effect. A new bump 
occurs each time the interpolator in the horizon finder changes the set of grid zones it uses. 
It is only visible because of the small range of Mah /Madm which we are viewing. 

In Fig. we show a similar graph for run k2, corresponding to a Kerr black hole with 
rotation parameter a/m = .68. Despite the higher rotation, all features noted for run kl 
are preserved to a high degree. Because there is more angular variation in the metric there 
is slightly less agreement between Mah and Madm- However, the pure Kerr spacetimes do 
pose more numerical problems than the distorted spacetimes, as noted in Paper I. In run k2 
the axis instability develops sooner at the highest resolution, and the results are unreliable 
after about t = 50M. 



C. Distorted Bowen and York 

In this section we discuss results for the six initial data sets labeled rO - r5. This repre- 
sents a sequence of distorted black holes with increasing values of the rotation parameter, 
a/m, beginning with a nonrotating case, labeled rO. 

The first case we discuss is run rO, which has the same construction as the "Brill wave 
plus Black Hole Spacetime" discussed in detail in Ref. This spacetime was evolved 
with an antisymmetric lapse across the black hole throat, so although the data set is in the 
same class as those evolved in Ref. it has been evolved with a new code capable of 
handling rotating black holes, and with a different slicing condition. The same simulation 
was discussed in Paper I, where a comparison of metric functions was made to the evolution 



obtained with the code described in Ref In Fig. |15]and Fig. [T6| we show the now familiar 
C. = 2 and £ = 4 waveforms extracted from the evolution, with a fit of the two lowest 
quasinormal modes in each case. The match is excellent. It is also important to note that 
no odd-parity radiation is present in this system and there are no £ = 3 or £ = 5 waveforms 
to show. ( It is important to point out that odd-parity does not mean odd-£, and even- 
parity does not mean even-£. It simply happens that with equatorial plane symmetry, which 
is present here, there are no odd-£, even parity modes, nor are there even-£, odd-parity 
modes.) The new dynamical variables used in this evolution remain exactly at zero when 
there is neither odd-parity distortion nor rotation. 



In Fig. O we show the mass of the apparent horizon for run rO, as defined in section [III A 



above. Note that initially it does not change, as the antisymmetric lapse freezes the evolution 
at the throat where the apparent horizon is found on the initial data set. Then at about time 
t = lOM, the horizon is found out at a larger radius and evolves slowly outward, its mass 
increases as more gravitational wave energy crosses the horizon surface into the hole. (This 
"jumping-out" of the horizon is a common property of all simulations we have performed 
with an anti-symmetric slicing condition as discussed in section |II1 A| above.) After a time 



of about t = 25M, there is a slow upward drift in the mass due to difficulties in resolving the 
peak in the radial metric function A that develops near the horizon . By this time, nearly 
all of the in-going gravitational wave energy has entered the black hole, and we compute the 
final mass of the hole to be Mah = -^^^Madm- We also have computed the total energy 
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radiated through the dominant i = 2 mode, as shown in Table |^ and plot it in Fig. |T7. 
Note that all the energy is accounted for to a high degree of accuracy. The total energy 
radiated by the black hole, as computed by integrating the Zerilli function, plus the final 
mass of the black hole, as determined by the mass of the apparent horizon, gives the total 
ADM mass of the spacetime, as one would expect. 

Next we turn to the run labeled rl. This calculation is similar to the previous run except 
that a noticeable amount of energy is now being radiated in the £ = 3 quasinormal mode, as 
shown in Fig. |18[ The value of a/m is extracted from C^, demonstrating that this technique 
is effective for small values of a/m, as shown in Fig. |T^ (although a/m = .35 is considered 
to be "small" in this study, the vast majority of black holes in the real universe are probably 
smaller than this [p^ ). 

Both runs r2 and r3 were typical of many simulations. Each provides waveforms that 
fit equally as well as r^, each locks on to the appropriate limit slice, and in both cases the 
radiation energy plus the final apparent horizon mass add up to the ADM mass. We show 
the plots for the radiation in the £ = 2, 3, 4, and 5 modes for run r2 in figures Figs. [2^, 
pT| , [2^ , and |2^ respectively. In each case we fit the wave function to the fundamental and 
first harmonic mode for the appropriate value of I. The value oi a/m used in the frequency 
fits was derived from the value extracted from the apparent horizon, although it made very 
little difference whether the dependence of the QNM frequencies on the rotation parameter 
was accounted for. Except in the case of extreme rotations, the QNM frequencies depended 
only weakly on a/m [^,^. In the case of the even parity radiation we include a constant 
offset in the fit to account for the non-sphericity of the spacetime. The fit is good in all 
cases. 

The calculation r4 is the second fastest rotator in this series. This rotation parameter 
{a/m = .7) was attained by reducing the strength of the Brill wave rather than by increasing 
the value of the angular momentum, and it is thus a different construction from runs rO to 
r3. Large positive amplitude Brill waves make it difficult to create spacetimes with large 
values of a/m, as they increase the mass significantly. We show the plots for the radiation 



in the i = 2, 3, 4, and 5 modes for this run in Figs. |2J, |2^, |2^, and In each case we fit 
the wave function to the fundamental and first harmonic mode for the appropriate value of 
i. Again, the fit is good in all cases. 

The run labeled r5 is a special case. No "Brill wave" distortion is present in the spacetime 
except the distortion necessary to make the hole conformally flat. Both the 3-metric and 
the extrinsic curvature are given by the simple solution of Bowen and York: He = 3J\Ef^ 
and Hp = 0. Bowen and York black holes are tricky to evolve, because they are so close 
to the Kerr black hole, yet this is the only method by which we have obtained a hole that 
rotates so rapidly. The presence of a Brill wave both stabilizes the spacetime and reduces 
the rotation parameter. 

We note that in all these spacetimes, as the rotation is increased, an ever smaller fraction 
of the energy available as radiation is able to escape from the black hole. The dominant 
energy source in this problem is the imposed gravitational wave at r/ = 1.0. Because of our 
coordinate transformation a constant rj value maps to an ever decreasing value of r as a/m 
approaches unity, and the width of the packet in the t] direction is effectively decreasing as 
a/m is increased, therefore less energy is distributed in the region outside the horizon. 

We note that in all these spacetimes it is approximately true that M^ad + MAHjinai is 
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equal to Madm- The "final" apparent horizon mass is, however, really the mass as given 
around 30 — 40M, since after that time large gradients in the radial metric variable A make it 
difficult to accurately obtain the horizon and bring about the well-known artificial inflation 
of the hole's mass at late times. 

VI. CONCLUSIONS AND FUTURE WORK 

Rotating vacuum black hole spacetimes are a significantly more complex type of ax- 
isymmetric vacuum spacetime than their non-rotating counterparts. In this paper we have 
applied a new code, designed to evolve rotating black holes, to evolve a new family of 
rotating black hole initial data sets. We developed a series of tools to analyze these space- 
times, including even- and odd-parity waveform extraction and various techniques to study 
the apparent horizons of these spacetimes. We showed that measurements of the horizon 
can be used to determine the mass, angular momentum, and oscillation frequency of the 
black holes. Distorted rotating holes are shown to oscillate about their stationary "Kerr" 
equilibrium configurations. 

Studies of a new class of odd-parity distorted non-rotating black holes were also made, 
and the waveforms and horizons were analyzed in detail. In all cases the normal modes of 
the black holes were shown to be excited and dominate the wave forms. We also determined, 
for the first time, the i = 5 quasinormal modes of a black hole by direct measurement of our 
numerically evolved spacetimes which were later verified by black hole perturbation theory 
techniques by Leaver. Finally, we were able to determine the final mass of the black hole 
by measuring the horizon area and compared it to the total ADM mass of the spacetime. 
The difference between these measurements agreed with the total energy radiated to within 
a few percent, showing the high degree of accuracy we are able to obtain in these studies. 

There is more work to be done in this for which our code will be useful. We plan to 
discuss, in future papers, other features of the initial data sets and to look at new data sets. 
We plan to study spacetimes with "cosmic-screw" type symmetry, which allows for rotation 
with equal and opposite amounts above and below the equatorial plane. Ultimately, we plan 
to study the collision of rotating black holes. 

More information about our group and the research we do is available through a world 
wide web server maintained by our group. We plan to have movies of some of the sim- 
ulations discussed in this paper available there. The URL for our group is fittp:/ /jean- 
luc.ncsa.uiuc.edu. 
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FIGURES 



FIG. 1. In this figure we plot the function Cr (the ratio of polar to equatorial circumference 

of the apparent horizon) as computed by our code for run labeled r4, a distorted rotating black 
hole. We have removed the early part of this plot where the horizon "jumps out" so that we can 
more clearly see the surface oscillations. The line labeled "equilibrium" is the offset from sphericity 
determined by our fit. This offset implies the system is oscillating about an "equilibrium" Kerr 
black hole with a rotation parameter of a/m = .70 black hole. 

FIG. 2. This figure shows Cr (the ratio of polar to equatorial circumference of the apparent 
horizon) as a solid line and an £ = 2 quasinormal mode fit form run r5. The horizontal short and 
long dashed line is the offset from sphericity, determined by the fit, showing that the black hole is 
oscillating about a Kerr hole with rotation parameter a/m = 0.87. 

FIG. 3. This figure shows the 3-space embedding of the distorted rotating black hole run 
labeled r4 at various times. Although the black hole's apparent horizon is initially prolate, it 
eventually settles down to the appropriate oblate Kerr black hole shape denoted by a dark solid 
line. 

FIG. 4. This plot shows the position at which the Euclidean 3-space embedding of the apparent 
horizon vanishes for run labeled r5, a Bowen and York rotating black hole, with angular momentum 
parameter J = 15 (solid line). As discussed in the text, for rapidly rotating holes the horizon 
embedding fails near the axis at an angle determined by the rotation parameter. The dashed line 
shows the angle at which the embedding fails for an "equilibrium" Kerr black hole with a/m = 0.87. 

FIG. 5. Wc show the horizon history diagram for the run labeled ol. The Gaussian curvature 
K is mapped to a grayscale on the surface of the horizon. The vertical axis gives the angular 
location on the horizon, and the horizontal axis traces out the time development. As discussed in 
the text, we see a pattern that repeats at twice the £ = 3 QNM frequency in accord with theoretical 
expectations. In the figure, dark regions are more highly curved than light regions. 

FIG. 6. In this plot we show the horizon history diagram for the run labeled o2. Hence the 
Gaussian curvature shows a diamond pattern, a result of the mixing of several frequencies of the 
£ = 3 and i = 5 quasinormal modes. In the figure, dark regions are more highly curved than light 
regions. Sec text for details. 

FIG. 7. We plot the extraction of the £ = 3 waveform from the odd-parity distorted non-rotating 
hole, run labeled ol. The dashed line shows the fit of this waveform to the two lowest £ = 3 
quasinormal modes. Practically all of the radiated energy is in this mode. 

FIG. 8. This figure shows the numerically extracted £ = 5 waveform and its fit to the two 
lowest £ = 5 quasinormal modes for the run labeled ol. 
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FIG. 9. In this figure we compare the £ = 3 wave form extracted from the (f) shift /3"^ (dashed 
Hne) and from the metric variable F = ge(f>/^'^ (sohd hne), where ^ is the conformal factor. 
Although we generally use the extraction based on the 3-metric, this demonstrates that other 
techniques can be used. Note that the waveform extracted from the (3'^ was normalized to have 
the same lower bound as the shift extracted from F. 

FIG. 10. We show the extraction of the £ = 5 waveform from an odd-parity distorted 
non-rotating hole, run labeled o2. It is clear that the fit (dashed line) agrees well with the numer- 
ically extracted function (solid line). 

FIG. 11. This figure shows a Fourier transform of the data extracted from the run labeled o2. 
The dotted line shows real part of the I = b frequency, computed via black hole perturbation 
theory. 

FIG. 12. This figure shows the apparent horizon mass defined in the text, for the run labeled 

ol normalized by the ADM mass (solid line). After about t = 40M the area of the horizon slowly 
drifts up due to difficulties in resolving metric functions, causing an over-estimation of the horizon 
mass. (See text.) The short-dashed line is drawn at 1 representing the total mass of the spacetime. 
The gap between the short-dashed and long-dashed lines is the amount of energy emitted through 
radiation. As one can see, the mass of the horizon and the energy of the radiation add up to the 
ADM mass. 

FIG. 13. This figure shows the apparent horizon mass for run kl [a/m = 0.48) at three reso- 
lutions. It is clear that the apparent horizon mass becomes more constant and agrees better with 
the ADM mass as resolution is increased. At all but the lowest resolution we have about 0.5% 
accuracy in the apparent horizon mass at the late times shown. 

FIG. 14. This figure shows the apparent horizon mass for run k2 {a/m = 0.68) at three reso- 
lutions. It is clear that the apparent horizon mass becomes more constant and agrees better with 
the ADM mass as resolution is increased. At all but the lowest resolution we have about 0.1% 
accuracy in the apparent horizon mass. 

FIG. 15. This figure shows the numerically extracted £ = 2 waveform (solid line) and the least 
squares fit to the two lowest i = 2 quasinormal mode for this distorted non-rotating black hole. 
This run was labeled rO. 

FIG. 16. This figure shows the numerically extracted ^ = 4 waveform and its fit to the two 
lowest 1 = 4: quasinormal modes for the run labeled rO. 
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FIG. 17. This figure shows the apparent horizon mass defined in the text, for the run labeled 
rO normalized by the ADM mass (solid line). The short-dashed line is drawn at 1 representing 
the total mass of the spacetime. The gap between the short-dashed and long-dashed lines is the 
quantity of radiation emitted. The mass of the horizon and the energy of the radiation add up to 
the ADM mass. 

FIG. 18. This figure shows the numerically extracted i = 3 waveform and its fit to the two 
lowest £ = 3 quasinormal modes for the run labeled rl. 

FIG. 19. This figure shows Cr (the ratio of polar to equatorial circumference of the apparent 
horizon) as a solid line and an £ = 2 quasinormal mode fit from run rl. The long and short dashed 
line is the offset from sphericity, determined by the fit, showing that the black hole is oscillating 
about a Kerr hole with rotation parameter a/m = .35. 

FIG. 20. This figure shows the numerically extracted i = 2 waveform and its fit to the two 
lowest i = 2 quasinormal modes for the run labeled r2. 

FIG. 21. This figure shows the numerically extracted i = 3 waveform and its fit to the two 
lowest i = 3 quasinormal modes for the run labeled r2. 

FIG. 22. This figure shows the numerically extracted £ = 4 waveform and its fit to the two 
lowest £ = A quasinormal modes for the run labeled r2. 

FIG. 23. This figure shows the numerically extracted £ = 5 waveform and its fit to the two 
lowest £ = 5 quasinormal modes for the run labeled r2. 

FIG. 24. This figure shows the numerically extracted £ = 2 waveform and its fit to the two 
lowest £ = 2 quasinormal modes for the run labeled r4- 

FIG. 25. This figure shows the numerically extracted £ = 3 waveform and its fit to the two 
lowest £ = 3 quasinormal modes for the run labeled r4- 

FIG. 26. This figure shows the numerically extracted £ = 4 waveform and its fit to the two 
lowest £ = A quasinormal modes for the run labeled r4- 

FIG. 27. This figure shows the numerically extracted £ = 5 waveform and its fit to the two 
lowest £ = 5 quasinormal modes for the run labeled r4- 
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TABLES 



Case 


J 


a/m 


kl 


2.5 


.481 


k2 


5.0 


.677 



TABLE I. This table gives summary data for evolutions of pure Kerr spacetimes. J is the 
total angular momentum of the spacetime, a/m is the usual rotation parameter in the Kerr metric. 



Case 


J 


Qo 


a 




n 


n' 


rO 


0.0 


1.0 


1.0 


1.0 


2 


N/A 


rl 


2.5 


1.0 


1.0 


1.0 


2 


N/A 


r2 


5.0 


1.0 


1.0 


1.0 


2 


N/A 


r3 


10.0 


1.0 


1.0 


1.0 


2 


N/A 


r4 


10.0 


0.5 


1.0 


0.0 


2 


N/A 


r5 


15.0 


0.0 


N/A 


N/A 


N/A 


N/A 


ol 


0.0 


2.0 


1.0 


2.0 


N/A 


3 


o2 


0.0 


2.0 


1.0 


2.0 


N/A 


5 



TABLE II. This table gives input parameters for each of the runs. J describes the amount of 
angular momentum in the system, Qq the amplitude of the Brill wave placed in the spacetime, a 
gives the width of the Brill wave, gives rjo the location of the peak of the Brill wave distortion, 
and the parameters n {n') describe the radial dependence of the distortion in the metric (extrinsic 
curvature) . 



Case 


{a/m)min 


{a/m)max 


{a/m)cr 


rO 


.000 


.000 


.000 


rl 


.342 


.600 


.351 


r2 


.432 


.698 


.436 


r3 


.511 


.759 


.512 


r4 


.696 


.827 


.703 


r5 


.868 


.886 


.870 


kl 


.480 


.481 


.481 


k2 


.675 


.677 


.677 



TABLE III. This is a table of a/m values calculated from the initial data. 

{a/m)min = J/^ADM initial slice and will be the final value if all the energy in the 

spacetime goes into the horizon, and {a/m)max = J/^^ah initial slice which will be the 

final value if the surface area of the horizon does not increase during the evolution. The parameter 
{a/m)cr is the angular momentum extracted from the horizon. 
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Case 


M^=2 
rad 


M^=3 

raa 


raa 


M^=5 
raa 


rad 1 J\UNi 


rO 


94.0% 


0.0% 


6.0% 


0.0% 


8.37e-2 


rl 


94.9% 


2.0% 


2.8% 


.3% 


3.40e-2 


r2 


94.8% 


3.4% 


1.5% 


.3% 


1.84e-2 


r3 


93.9% 


5.2% 


.7% 


.2% 


9.41e-3 


r4 


92.1% 


7.0% 


.6% 


.2% 


8.30e-3 


r5 


95.7% 


4.1% 


.2% 


.02% 


l.OOe-3 


ol 


.05% 


99.9% 


.02% 


0.0001% 


1.06C-4 


o2 


.06% 


78.9% 


0.00% 


21.0% 


2.88e-4 


kl 


N/A 


N/A 


N/A 


N/A 


4.62e-8 


k2 


N/A 


N/A 


N/A 


N/A 


2.00e-7 



TABLE IV. This table shows the total energy emitted as radiation divided by the ADM mass in 
the right-most column, and it shows the fraction of this energy emitted in the first four quasinormal 
modes of the black hole spacetime. 
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http://arXiv.org/ps/gr-qc/9412073vl 



1=2 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=2 




-2.0 ' ' ^ ' ^ ' ^ ' ^ ' ' 

0.0 20.0 40.0 60.0 80.0 100.0 

Time (M^,J 



This figure "fig2-15.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9412073vl 



1=4 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=4 




Time (M^,J 



This figure "fig2-16.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9412073vl 



Apparent Horizon Mass vs. Time 




This figure "fig2-17.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9412073vl 



/=3 Wave Mode Extraction 



Fit to Fundamental and First Harmonic of 1=3 

0.20 I 




100.0 



Time (M^,J 



Oscillating Horizon for run rl 




1=2 Wave Mode Extraction 



Fit to Fundamental and First Harmonic of 1=2 




0.0 20.0 40.0 60.0 80.0 

Time (M^,J 



0.20 



0.10 



0.00 



-0.10 



-0.20 



0.0 



/=3 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=3 




20.0 



40.0 60.0 
Time (M^,J 



80.0 



100 



1=4 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=4 



0.10 




1=5 Wave Mode Extraction 



Fit to Fundamental and First Harmonic of 1=5 

0.060 I 




-0.060 ' 

0.0 20.0 40.0 60.0 80.0 

Time (M^,J 



1=2 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=2 




Time (M^,J 



/=3 Wave Mode Extraction 



Fit to Fundamental and First Harmonic of 1=3 

0.20 I 




-0.20 ' ^ 

0.0 20.0 40.0 60.0 80.0 

Time (M^,J 



0.04 



0.02 



0.00 



-0.02 



-0.04 



0.0 



1=4 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=4 




20.0 



40.0 
Time (M^,J 



60.0 



1=5 Wave Mode Extraction 

Fit to Fundamental and First Harmonic of 1=5 




